On the ChristofFel-Darboux kernel for 
random Hermitian matrices with 
external source 

Jinho Baik * 

September 23, 2008 

Abstract 

Bleher and Kuijlaars, and Daems and Kuijlaars showed that the cor- 
relation functions of the eigenvalues of a random matrix from unitary 
ensemble with external source can be expressed in terms of the Christoffel- 
Darboux kernel for multiple orthogonal polynomials. We obtain a repre- 
sentation of this Christoffel-Darboux kernel in terms of the usual orthog- 
onal polynomials. 



1 Introduction 

Fix a Hermitian matrix A, and a function V{x) on M, which decays sufficiently 
fast. Consider the ensemble of n x n Hermitian matrices with the density func- 
tion defined by 

p(M)-c„-e-T'W*^)-^^^). (1) 
Here c„ is the normalization constant. 

c„= (^le-T't^t^^'-^^^^dA/^ (2) 

where the integral is over all n x n Hermitian matrices. Due to the unitary- 
invariance of the density function, we can assume without loss of generality that 
j4 is a diagonal matrix: 

A = diag(ai,a2, • • • , a„) (3) 

where ai > 02 > • ■ • > a„. When A ~ 0, the above ensemble is the standard 
unitary- invariant ensemble in the theory of random matrices |15j . The case 
when A 7^ is called a random matrix model with external source, and has 
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been studied in, for example, [TH [1^1 [21 [TJ [S]. It is well known that, using 
the famous Harish-Chandra-Itzykson-Zuber integral [131 HIIj © induces the 
probability density function of the eigenvalues Ai > • ■ • > A„ of M given by 



p(A) = C„ • det( 



l<i<j?'<n — 1 



(4) 



where C„ = C„(ai, • • • , a„) is the (different) normalization constant. 

In random matrix theory, one of the fundamental object of interest is the 
correlation functions of the eigenvalues. When ^ = 0, the correlation functions 
are expressible in terms of the Christoffel-Darboux kernel of the orthogonal 
polynomials with respect to the weight e"^^^) [TS]. Using the Ricmann-Hilbcrt 
problem (RHP) representation of orthogonal polynomials [H], the asymptotics 
of orthogonal polynomials with respect to a general class of weights have been 
obtained [9l llOj. hence yielding limit theorems for correlation functions. 

When A 7^ 0, it was shown by Bleher and Kuijlaars [1| that the correla- 
tion functions of the eigenvalues can be expressed in terms of a kernel involving 
a sum of multiple orthogonal polynomials. This summation was further sim- 
plified by Daems and Kuijlaars [7], yielding a Christoffel-Darboux formula for 
multiple orthogonal polynomials. The purpose of this paper is to re-express 
the Christoffel-Darboux kernel obtained by Daems and Kuijlaars in terms of 
orthogonal polynomials. 

In order to state the result, we first discuss the result of [7]. Assume that 
n~ r eigenvalues of A are zero. If A has no zero eigenvalue, then we can simply 
shift ^ by ^ — ail and V{x) by V{x) — aix. Assume further that all non-zero 
eigenvalues of A are simple: 



A = diag(ai,a2, ■ • • ,0^,0,0, ••• ,0) 
where ai > a2 > ■ ■ ■ > ar and aj ^ for all j = 1, • • • 



(5) 

Where there are 

multiple eigenvalues, one obtain obtain the relevant result by taking limits such 
as ai — > a2 in the algebraic formulas below. 

Let Z be the solution to the following {r + 2)x (r + 2) RHP: 

• Z(z) is analytic for z e C \ M and is continuous up to the boudary, 

• Z+{z) = Z_(z)Jz(z) for z e M where 



3z{z) 



/I 





,-V{z} -V(z)+axz 



yo 



g-y(z)+a,.z\ 






1 / 



(6) 
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Z(z) = {lr+2+ 0{Z-^)) 








z-^ ■ 












■ 


■ --7 



(7) 



Here for 2 G M, Z-|_(z) (resp. Z_(z)) is the limit of 7i{w) as w ^ z where w is 
in the upper (resp. lower) half plane. It is known that the solution Z exists, is 
unique, and is expressed in terms of multiple orthogonal polynomials of type II 
[TTI. Define the kernel 



-V{x) 



2Tri{x — y) 



1 



/1\ 





(8) 



Note that since (1,0, ■ • • ,0)* is a right eigenvector of Jz{x) corresponding to 
the eigenvalue 1, and (0, 1, e^^^, • ■ • ,e°'^y) is a left eigenvector of Jz{y) corre- 
sponding to the eigenvalue 1, the kernel IC{x,y) is independent of the choice of 
the plus or the minus limits of Z(.t) and Z(y). Daems and Kuijlaars [7] showed 
that the correlation functions of the eigenvalues of random Hermitian matrix 
with density ([T]) when A is of form ^ is expressible in terms of K,. For example, 
the largest eigenvalue has the distribution 



P(Amax <s)=dct(l-/C|(,,oo)) 
-2/ 



(9) 



where IC \(s,oo) is the operator on L {{s, 00)) whose kernel IC{x, y). 

Let us compare the kernel IC with the Christoffel-Darboux kernel for orthog- 
onal polynomials. Let ICq be JC when A = 0, i.e. without any external source. 
Concretely, /Co is expressed in terms of the following RHP. Let X be the solution 
to the following 2x2 RHP: 

• X(z) is analytic for z G C \ M and is continuous up to the boudary, 

• X+(z) = X_(z)Jx(z) for z e M where 



Jx(^) 



1 e-^(^) 
1 



(10) 



and 



as z ^ 00. 



X{z) = {h + 0{z-'))(l 



(11) 



^Throughout the paper I/; denotes the k X k identity matrix. 
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This is precisely the RHP for orthogonal polynomials [12] ■ Let 7r„(z) = z" + • • • 
be the monic orthogonal polynomial of degree n with respect to the weight 

^-V{x). 

nn{x)Tr„,{x)e-^^''Ux = 0, n ^ m. (12) 

Let Kn be the positive constant such that {KnT^niz)} becomes a sequence of 
orthonormal polynomials. Then |12| . 

^^"^^(-2^*4"it-iw *)' ^^^^ 

where the second column is given by the Cauchy transform of the first column 
times the weight e~^^'^\ However we do not need the explicit form of the second 
column throughout this paper. Theifl 

K^o[x,y) = 

""-y (14) 



i)x+(y)-ix+(x)(i 



2Tri(x — y) \0 

Now we state the result of this paper. Define the column vector of length r, 
t(z) = (7r„_r(2;) 7r„_,.+i(z) ••■ 7r„_i(z))*. (15) 
Also define the column vector 

v(z) = (e"i^ e"'"' ■■■ e""-^)* (16) 

and set 

/•OO 

w{zf^w{yY~ ICo{s,y)v{syds. (17) 



Define an r x r matrix 

/OO 
t(.s)v(s)*e~^('')ds. (18) 
-OO 

Theorem 1. Let A be as in ((5)). Then the matrix B is invertihle and 

(/C(x,y)-/Co(x,2/))e^(-) =w(y)*B-it(x). (19) 

From a standard linear algebra computation (see ([59]) below), the above 
result can also be written as 

(/C(x,y)-/Co(x,2;))e^(-' = -^dct^^ ^g)*) . (20) 



■^Usually /Co is defined by the term e ^(^) replaced by e 2 . But they are 

conjugate of each other, and have the same spectrum. 
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Note that the RHP © for multiple orthogonal polynomials are of size higher 
than 2, contrary to the RHP (fTO)) for orthogonal polynomials which is of size 
2. There have been many exciting recent progresses in solving k x k RHP 
asymptotically when fc = 3,4 (see for example, [21 II [51 [51 [TT] ) . Nevertheless the 
asymptotic study of RHP of general higher size is still a challenge. The main 
motivation of this paper is to find an expression of the correlation function 
involving only 2x2 RHP's. 

A special case is the so-called spiked model of random Hermitian matrices 
with external source when n tends to infinity but r is kept fixed. In this case, the 
formula (fTO)) seems well suited for the asymptotic analysis since the right-hand- 
side of is of fixed rank r. In order to control this term asymptotically, one 
needs the first r sub-leading terms of orthogonal polynomials in the asymptotic 
expansion explicitly. This requires a refinement of the analysis [H [10] . In the 
future work, we plan to work on this problem for general weight V{x). The 
particular case when V{x) = and also the Laguerre version of the problem 
have been studied in [51 [TB] in the context of sample covariance matrices, which 
obtained various limit theorems and interesting phase transition phenomenon. 
The goal of the future work is to prove the universality for spiked model with 
general weight V. On the other hand, when the rank r oi A grows with n, (|19p 
does not seem to be suitable for asymptotic analysis as the rank of the right- 
hand-side of (|19p grows with n. For such case, one needs to analyze of the 
RHP © for multiple orthogonal polynomials directly, and [21 [H [S] worked out 
such asymptotic analysis when V{x) = x^, and the half of the eigenvalues of A 
are zero and the other half are the same, say a. 

The proof of the theorem is based on algebraic manipulations of RHP's. 

Acknowledgments: The work of the authors was supported in part by NSF grants 
DMS 0457335, DMS 075709. 

2 Proof of Theorem [1] 

Set 

M(z) := Z(z) 0). (21) 

Then M(z) solves a new RHP. Let denote the (i, j)-entry of X and set 

-(-)=fx11) = f 9 -"^'^ ( (22) 
\X2i{z)J y-27riK^_i7r„_i(z)y 

The matrix M(z) satisfies: 

• M(z) is analytic for z G C \ M and is continuous up to the boudary. 
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M+(z) = M_(z)Jm(^) for z e M, where 



[\ Xii(2)e-^(^)+'^i^ 
1 X2i(z)e-^(^)+'^i^ 
1 



Vo 

I2 x(z)v(z)*e-^(^) 
L 



Xii(z)e-'^W+'^'-^\ 
X2i(z)e-^(^)+''- 




(23) 



and 



00, 



M(z) = (l,+2 + 0(z-i)) 



/ 


Z"X22(Z) 


-^"Xi2(z) 


■ 


■ \ 




-2-"+'^X2l(z) 


z-"+'-Xn(z) 


■ 


















V 








• 





(24) 



Here 1^+2 + 0(z~^) in (|24p and 1^+2 + 0(z~^) in ([7]) are the same function. It 
will be clear in the below that the second column of this term plays a distinct 
role, and we denote it by 



(0,1,0,--- ,0)* + (ei(z),e2(^),<5i(;2),--- ,-5^(2))* 
O(z-i) equals Z(z) diag(z", z"!, 



(25) 



Since the term 1^+2 + 0{z^^) equals Z(z) diag(2;", 2""+'', ■ • ■ ^z~^)~^ and 
Z(z) is given by the multiple orthogonal polynomials of type II [17], the term 
0(z~^) indeed have an asymptotic scries in the powers of . Let 



e,(z) = eiz-' + 4^-' + • • • + ,lz-^ + 0{z-^-'), 
4(2) = 5\z-^ + ^Iz-^ + • • • + ^Iz-'- + 0(2-'-i), 



fc= 1,2, 
fc = 1, • • • 



(26) 



as 2 ^ 00 for some constants e j. , j ~ \ 



fc, J - • • • = 1,2 and £ = 1, - • ■ 
We solve M in terms of X. Write M in a block-form: let 



M(z) 



M(ii)(2) M(i2)(2; 

M(21)(2) M(22)(z; 



(27) 



where M^^-^^ denotes the upper-left 2x2 block of M, and M^^^^ denotes the 
lower-right r x r block of M. The jump condition on M then becomes 



m%^'^{z) yi%^\z)\ _ (m^^^\z) ML"^(z)x(z)v(z)*e-^(^) -|-ML^^'(z)\ 
M%^\z) yi%^\z)) ^ \m^l^\z) ML^^^(z)x(z)v(z)*e-^(^) +ML^^'(z)y 

(28) 
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forz e R. Hence M(")(z) is continuous on R and therefore is entire. As z — > co, 
the RHP for X impUes that z"X22(z) = 1 + 0(z"^) and z"Xi2(z) = 0{z-^), 
and hence 

l + 0(z-i) ei(z) \/ z"X22(2) -2"Xi2(2) 
O(z-i) l + e2WM-2-"+''X2i(z) z-'^+'-XnW 



1 




+°<-')+(i:'ll))<- 



'^X2i(z), z-"+'-Xn(z)) . 

(29) 



Let and {cj} be the coefficients of Xn and X21: 

Xn(^) = 7r„(z) = 60^" + biz"-^ + b^z''-^ + ■ • • + 
X2i(2;) = -27riK^j_i7r„_i(z) = cqz" + ci; 
Note that 

60 = 1 Co = 0. 

From tlic Liouville's theorem, M^^^^z) is a polynomial. This polynomial is 
obtained by expanding out the expansions of ekiz) and Xij(z) in (j29p . and 
collecting only the terms of non- negative powers of z. Explicitly, set 

Pi{z) boz' + biz''^ + b2z'-^ + --- + b-j 



(30) 
(31) 



Cqz' 



- C\Z- 



J-1 



C2Z 



Then 



where 



M("'(z) =D(z) 



D(z) := 



1 E(z) 



1 



and 



(-Q\-\z) Pl-Hz)\ 

-Ql-\z) p:^-Hz) 

E(z):= . . 

V -Quiz) P^{Z) ) 
The jump condition for M^^^) (ggg 

]V[^^^^(z) mL^^'(z) + M(")(z)x(z)v(z)*e-^("\ z e 



(32) 

(33) 
(34) 

(35) 

(36) 



Here mL"^ = M(") since M(^^'(z) is entire. Using the Plemelj's formula for 
the solution of an additive RHP, we obtain 



M^^^\z) = — [ ^M(")(s)x(s)v(s)*e-^(^)ds, 
27ri Jts,s-z 



(37) 
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for z G C \ R. Note that there is no additive term of an entire function in this 
formula due to the asymptotics 



M(12)(z) = 0(z-2)^ 



(38) 



which follows from ([24)1 . This asymptotic condition, applied to ([37|) . further 
implies that 



M(")(s)x(s)v(s)*e-^(^)ds = 0. 



From ([35|) . the equations become 



"2 "=2 



1 Q = C, 



where 



and 



Q := / E(s)x(s)v(s)*e-^(")rfs. 
Jr 



(39) 

(40) 

(41) 
(42) 



D(s)x(s)v(s)*e-^('')ds. 

We will sec ((50)) below that Q is invertible. Hence (|43)) and p3|) imply that 

M(")(2) = D(z) + CQ-1E(z). (43) 

Now consider M^^^^. The jump condition ([28]) implies that M(^^)(z) is an 
entire function. As z — * oo, from (|24p . we have 



M(2i)(z) 



/O(z-i) <5i(z)\ 



\Sr{z)J 



Z"X22(Z) -z"Xi2(z) 

-z-"+'-X2i(z) z-«+'^Xn(z; 



(-z-"+'^X2i(z), z-"+^Xniz)) 



(44) 



From the Liouville's theorem, M^^^^z) is a polynomial of degree at most r — 1, 
and is given by (recall (US])) 



M(2i)(z) = 





■ 5\\ 




■ 51 


VI 51 ■ 


■ 5lj 



E(z), 



(45) 
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where £(2;) is defined in ([55]) . 

The jump condition for M^^^^ is 

M^^^' {z) = mL^^' (z) + M(21) (z)x(z)v(z)*e-'^(^) , 2 G I 
Using the asymptotic condition 

M(22)(2) = ^-il,, + 0(z-2)^ z^oo, 
the Plemelj's formula implies that 

M(22)(z) = / ^M(2i)(s)x(s)v(,s)*e-^(^)ds, z e C\: 
The asymptotic condition (l46l) further implies that 

/ M(2i)(s)x(s)v(s)*e-^(^)cis = -2rf,. 
From (pS]) . this implies that 



(5\ 51 



5l\ 
S'2 



Q = -27ril^ 



where Q is in ((4T|) . This implies that Q is invertible, and 

M(21)(z) = -27riQ-iE(z). 

Therefore, combining dST]), gS]), ([48]) and ([51]), the matrix M(z) is 
given by 

_ (Biz) + CQ-iE(z) (CiD)(z) + CQ-i(CiE)(z) 
ivil^j - I _2^iQ-iE(z) -27riQ-i(CiE)(z) 



I2 C 

-27ra, 



\ /D(z) (CiD)(z) 



Q-V (CiE)(z) 



where we use the notation 



2711 s - z 

for 2 X 2 matrix-valued function F{z) on M. 



Now consider the kernel ([5]). From (HI]), 
27ri/C(x,y)e^(^)(a;-y) 



(0 1 



x+(j/)-^ 





X+{x) 

Ir 







On the other hand, from p^ . 

27ri /Co(x,?;)e^(^)(x-y) 
= (-X.,(,) Xn(.))fJ;;g 



/Xii(a;)\ 
X2i(a;) 




V / 



(54) 



= (-X2i(y) Xii(?/) 

Hence 

27ri(/C(a;,y)-/Co(a;,y))e^(^' 



aar-y) 



/Xii(x)\ 

X2i(.t) 



V ; 



(55) 



x-y 



-X2i(y) Xniy) e-^y ••• e'^-'') (M+(y)-iM+(x) 



/Xn(x)\ 
X2i(a;) 




V J 



(-X2i(y) Xn(y) e''^^' ••• e-^y)M+iy) 



_iM(x)-M(y) 



X — y 



/Xii(x)\ 
X2i(a;) 


V ; 



(56) 

Here in the last line, we replaced ^+'-^'>~^+^y'> by since only the 

' ^ x—y ^ ^ — y 

first two columns arc needed when multiplied by the last vector, and the first 
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two columns of M(z) arc entire. Using the formula ([5^ for M, we obtain 
27ri(/C(a;,y)-/Co(a;,y))e^(-) 

^D(y) iClB)iy)Y' ^""^^^^ 



-X2i(y) Xii(y) e° 



E(y) (Ci.E)(y); I e(.)-e(,) 



x(a;). 



(57) 



Here note that since M(z) is invertible and det M(z) = 1, the last matrix in (|52p 
is invertible and 



det 



D(z) (CiD)(z) 
E(z) (CiE)(z) 

Using the general identity 

t A-i 1 



■ det 



-27ri)-''detQ. 



w* 

W2 A 



det A 

for an invertible matrix A and vectors wi,W2, we obtain 



27rz(/C(a;,y)-/Co(a;,y))e^(^) 



{-27riY 
detQ 

{-2niY 
detQ 



det 



det 





E(:rO-Efa) 



x(y) v(2/)* 

^(^^x(x) D(2/) (CiD)(y) 
\^M£i^x(z) E(y) (CiE)(y) 



X2i(j/) Xn(y) e-i^' •• 

D(y) (CiD)(y)\ 
E(y) (C|E)(y); 



where 



and 



-1 

1 



x(y) = 



v(z) 



□aiz a2Z 



-X2i(y) 
Xii(2;) 



as defined in pB)) . Set 



A := det 



x(2/) v(2/)* 

5M£)x(a;) E(y) {CiE){y) ^ 



so that 



2TTi{lC{x,y) - /Co(.T,?;))e 



y(2;) _ 



(-27ri)'' 
detQ 



A. 



(58) 



(59) 



(60) 
(61) 

(62) 

(63) 
(64) 
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We have 










x(y) 


v(y)* 


A = dct 


D(x) 


^x(x) 


D(y) 


(C|D)(y) 






:f'^)x(.) 


E(y) 


(C|E)(y) 




/ 







v(y)* - 


= det 


D(x) 

X 




D(y) 


(CiD)(y) 




Eix) 
\ X 


EfMx(a;) 


E(y) 


(CiE)(y) 




(65) 



(13)-block is 



2-Ki{s—y) 



x(y)* 



1 



27ri 7r (s - y)+ 
K-o{s,y)w{sfds = w(y)* 



(s)v(.s)*e-^('*)ds 



v(y)*-x(y)(C|l2)(y) = v(y)* 

(see (HZ])). On the other hand, since the first row of D(z) is (1, 0) and 

D(,s) - D(y) 



(66) 



(CiD)(2/)-D(j/)(C|l2)(y) = -^ / 



-x(s)v(.s)*e-^(*)ds, (67) 



the second row of the matrix in (p5|) equals (0, 1, 0, 0, • • • ,0). Consider the last 
row of the matrix. Since the last row of E(z) is (— Q° (z), P°(z)) = (— cq, feo) = 
(0, 1), we find that the last row of the matrix in ((65)) is (0, 0, 1, 0, • • • ,0). Hence 
the determinant of the matrix in (j65|) is unchanged, up to the factor (—1)'^, if 
we remove the second and the last row (which is the (r + 3)rd row), and the 
second and the third column. Therefore, noting that the second row of D(z) is 
same as the first row of E(z) with r — 1 replaced by r, we arrive at the formula 

/ vf{yf \ 

^=^-^y''-\ii^^i.) (ciF)(,)-F(,)(c|i,)(,)l 



where 



F(z) 



/ ~Ql{z) P^z) \ 



(I z . 
1 • 



\0 



\ / -Cr br \ 

-Cr-1 br-1 



J \ -Co bo J 



(69) 
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The last line is an r x (r + 1) matrix times an (r + 1) x 2 matrix. Now 
F{s)~F{y) 



s-y 



/O 1 s + y 
1 



yo 



s''-l + s'-2y + .. 



/ -Cr br \ 

-Cr-l br-1 
-Cr-2 br-2 



S + y • • 



/I 

1 
\0 

U(y)U(s) 



V -Co bo ) 



-Cr-2 br-2 



where 



U(z) := 



= U(2;)E(s) 
/I z ... 

1 ... 

\Q ... 



-,r-l 



1 / 



Therefore we obtain 



A = (-l)''det 

= (-!)'■ det 

= (-l)Met 

since det U(y) = 1. 
Now we compute 







(u(2/)E(x)x(x) 




uM 



w(y)* 

/«E(.)x(,s)v(5)*e-^(^)d. 
w(y)* 



E(x)x(z) ^ /j, E(s)x(s)v(s)*e-^(^)ds 

w(y)* 
E(x)x(x-) 



E(z)x(z) 



/-gri(z)Xn(z) + W2i(^)\ 

-g;;-2(z)Xn(z) + pr'WX2i(^) 



-gi(z)Xn(z)+P„i(z)X2i(z) 
V -gO(z)XnW+P°(z)X2i(z) ) 
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The polynomial 

-Oi(z)XnW + P^(z)X2i(z) 

= -{caz^ + ■■■ + Cj){boz'' + ■■■ + 6„) + {bozJ + ■■■ + 6j)(coz" + ••• + €„) 
= -{cobj+i - boCj+i)z"^^ H 

(74) 

is of degree n — 1. Hence it is a linear combination of the orthogonal polynomials 
TTo{z), • • ■ , 7r„_i(z). Define the real inner product, 



(75) 



Recalling that Xii(z) = 7r„(2:) and X2i(z) = — 27riK^_j^7r„_i(z), we have for 
fc = 0, 1,-- - ,n- j-2, 

< -Q^Xii +P^X2i,7rfc >= - < TTn^QiTTk > ~2mn'^^_^ < 7r„_i,P^7rfe >= 

(76) 

since Q^tt^ is of degree k + j — 1 (recall that cq — 0) and P-^TTk is of degree k + j. 
Also since 6o = 1, 

< -Q^^Xii + P,-JX2i,7r„_j_i > = -27riK^_i < 7r„_i, P^7r„_j_i > 

= -2niKl_^ < 7r„_i,z"~^ > 

= -27rz. 



(77) 



Noting that < Tr^jTrfe >= Kj,^, we obtain 



E(z)x(z) = -27riSt(z) 



where 





V 



(78) 



(79) 



<-i/ 



A = (-l)''det 
= det S det 



and t(z) is defined in (fT5|) . From (|41]l . this implies that 

Q = 27riSB 

where B is defined in (fT8|) . As Q is invertiblc, so is B. Thus 

w(y)*' 
-27riSt(a;) -SB^ 

w{yY' 
^2'Kit{x) B 

using dSSl). As detQ = (-2771)*^ det(SB), we obtain from §^ that 

2TTi{lC{x, y) - ICo{x, y))e^(''' = 27ri\v{yYB~h{x). 
This completes the proof of the Theorem. 



= -27ridet(SB)w(y)*B~H(a;) 



(80) 



(81) 



(82) 
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